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A novel mechanism is proposed for single and double qubit 
state manipulations in quantum computation with four-fold 
degenerate energy levels. The principle is based on starting 
with a four-fold degeneracy, lifting it stepwise adiabatically by 
a set of control parameters and performing the quantum gate 
operations on non-degenerate states. A particular realization 
of the proposed mechanism is suggested by using inductively 
coupled rf-squid loops in the recently observed macroscopic 
quantum tunneling regime where the energy eigen levels are 
directly connected with the measurable flux states. The one 
qubit and two qubit controlled operations are demonstrated 
explicitly. The appearance of the flux states also allows pre- 
cise read-in and read-out operations by the measurement of 
the flux. 

PACS numbers: 



Recent advances on the experimental efforts to demonstrate 
the fundamental single and two bit quantum gates by ioa 
traptJ and nuclear magnetic resonance (NMR) experimentso 
as well as the observation of the coherent Rabi oscillations 
in Cooper pair boxesfj are hinting at new challenges await- 
ing the realization of the quantum computational devices. It 
has been shown theoretically that the basic gate operations 
based on the unitary transformations of the qubits. can be per- 
formed by single and double qubit manipulations!] . Therefore 
it is a necessary first step that a proposed mechanism should 
demonstrate these fundamental operations as well as a suc- 
cessful read-in and read-out before any decoherence comes 
into play. A next step which is not detailed in this letter is to 
construct arrays of independent qubits or coupled qubit pairs 
to perform parallel gate operations. In our suggested mecha- 
nism the state space is four dimensional corresponding to the 
state space of a pair of qubits. The paired qubit states are 
inseparable; nevertheless, each qubit can be manipulated in- 
dividually. The mechanism also allows controlled operations 
on selected bits rather easily without affecting those bits that 
need to be unchanged. A more detailed discussion of the pro- 
posed mechanism with its specific experimental realizations 
using rf-squid loops will be published elsewhere. 



A. The mechanism 

The mechanism is based upon coupling a pair of ideal two 
level systems with manifestly degenerate energy levels at the 
energy Eo by two complex coupling parameters as indicated 
in Fig. la. The coupling Ci lifts the four- fold degeneracy to 
two- fold (Fig. lb). The second coupling C2 couples the two 
remaining degenerate levels in each pair and lifts the de- 
generacy completely (Fig.lc). The states in the initial four- 
fold degenerate configuration at Ci = C2 — are labelled 
by \mn), (m,n = 0,1). We consider them to be the eigen 
states of some underlying non interacting system. In the 
qubit language below the first (second) bit is m(n). The 
states \mn) are the components of the four dimensional row 
vector (| 11), 1 00), |10), |01)) which is to be considered as the 



basis for the model Hamiltonian. The Hamiltonian includ- 
ing the degeneracy breaking couplings is given in this four 
dimensional basis by the Hermitian matrix 
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Without loss of generality we assumed that the couplings are 
real and positive and that they are applied adiabatically in 
the order C\ first and C2 second (the order is commutative) 
with < C2 < C\. The diagonalizing matrices for indepen- 
dent couplings are unitary transformations acting on the four 
dimensional basis. These are U\ for Ci 7^ 0,C2 = and U2 for 
Ci = 0, C2 7^ as given by 
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When both couplings are turned on, the Hamiltonian (|jj) is 
diagonalized by U —U2U1 asHj = diag{Eo+Ci+C2, E0+C1 — 
C2, E0—C1+C2, E0—C1—C2} which correspond to the levels in 
Fig.lc. For arbitrary couplings &,Ca 7^ the eigen energies 
are all nondegenerate and they are distributed symmetrically 
around the major center Eo and the two minor ones Eo ± 
Ci . The single and double qubit quantum gate operations are 
performed ideally by manipulating the strengths of at least 
one of the coupling constants and making use of the invariance 
of the major center under adiabatic changes of both couplings 
and of the two minor centers under that of C2. In Fig.lc we 
have indicated the eigenstates with up and down arrows. In 
the double-SQUID model considered below they correspond 
to the flux states. These states are connected with the initial 
basis states \mn) by 
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Note that the ordering of the \mn) basis is made with respect 
to the ordering of the energy eigen values of the flux states on 
the left hand side in which is different than the ordering 
used for constructing the Hamiltonian matrix (Q). For later 
convenience we introduce the short notation 
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B. An rf-SQUID realization of the mechanism 

Physically, this suggested four state mechanism can be re- 
alized by using inductively coupled rf-squids as depicted in 
Fig. 2. The SQUIDs must oncfiate in the macroscopic quan- 
tum tunneling (MQT) regimeQ'Q which requires temperatures 
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on the order of mili Kelvin for standard Josephson junctions. 
The four level system of Fig. la is constructed by preparing 
two identical SQUIDs A and B in which each qubit is rep- 
resented, to a good approximation, by its two lowest energy 
eigen states (denoted as |0) and |1}); in the MQT regime these 
are nearly degenerate (pseudo-degenerate) in energy (due to 
the fluctuations in the flux) compared to the next (third) level. 

A qubit which is close to being ideal is characterized by an 
isolated hypothetical two-level system with a high qubit qual- 
ity ratio rj = (E2—E1)/ (E1—E0), where Eo and E\ are the two 
lowest levels and E2 is the next (well-separated) level. In the 
case of manifest degeneracy, i.e. Eo = E\ 7^ E2 one obtains 
77 = 00, the ideal case. For the realistic case here using cou- 
pled SQUIDs, we search for a pseudo degenerate configuration 
(by performing a numerical search at different parameter val- 
ues) in which we obtain sufficiently high qubit quality factors 
on the order of r\ ~ 40. In Fig. 3 the calculated wavefunctions 
for the two lowest states of a symmetric double-well single- 
SQUID potential is plotted for the specific parameter values 
which are junction capacitance C = OApF, the loop inducn 
tance L — WOpH and the potential modulation parameters 
[3l — 2nI c I/&o =1.1 where I c is the critical Josephson cur- 
rent, I is the superconducting current in the loop and $0 is 
the superconducting flux quantum $0 = hc/2e. In the nu- 
merical calculations the range of the junction and the loop 
parameters are chosen such that the pseudo-degenerate eigen 
states are located in the double well, whereas the third level 
is above the central maximum of the double well potential. 
For the flux dominated regime the Josephson energy Ej is 
expected to be much larger than the charging energy E c in 
the junctions. Considering that the sum of these energies is 
fixed (Et ~ E c + Ej is the total energy in the noninteract- 
ing SQUID loop), we find for our parameters Ej/Et ~ 95% 
for the symmetric configuration of the double well potential. 
If the symmetric configuration is tilted, the additional level 
splitting raises the charging energy. For the range of interest 
of the bias flux (see Fig. 4) a drop to 90% in the above ratio 
is observed as a result of the tilting. The pseudo- degener- 
ate single SQUID states |0) and |1) have respectively even and 
odd parities. We will refer to these states as the one-bit parity 
basis. The odd and even superpositions of the parity states 
are given by | T > = (|0> + |1»/V2 and | J.) = (|1> - |0»/V2. 
These are the single-well localized flux states. Hence, we re- 
fer to them as the one-bit flux basis. The single-well local- 
ized means that these one-bit states actually refer to a single 
SQUID where they are well confined to a single well of the 
one-squid double-well potential. The calculated flux states 
corresponding to the chosen SQUID parameters above are 
also shown in Fig. 3 (the solid triangles). 

We now change from the single bit-single SQUID to the 
double bit- double SQUID configuration. The case for two 
identical SQUID loops in the pseudo-degenerate regime with 
no mutual interaction is given by Fig. la. Inductively cou- 
pling these two identical SQUID loops in Fig. 2 by a weak in- 
ductance M(<C L), lifts the four-fold degeneracy to two- fold 
degeneracy. In the coupled SQUID model the weak mutual 
inductance cannot be varied after fabrication. This mutual 
inductance corresponds to the coupling C\ (Fig. 1) and we be- 
lieve that M ~ 10 -2 L is a reasonable range to consider in the 
experiment. Once the SQUIDs are manifactured, M is fixed, 



therefore the four fold degenerate configuration in Fig. la is 
out of reach. However we still have leverage on the second 
coupling. This second coupling is the tilt of the symmetri- 
cal Josephson potential of the two SQUIDS simultaneously 
deriving both potentials to unsymmetric configurations by 
means of an additionally applied flux AQ ex (The total flux 
is $ C i = A& ex + $o)- This is modelled by the second cou- 
pling C2 which lifts the remaining two-fold degeneracy com- 
pletely. The combined effects of the mutual inductance, which 
we now consider to be M = 0.03 L, and the tilt can be cal- 
culated numerically by diagonalizing the inductively coupled 
two-SQUID Hamiltonian. The calculated eigenenergies are 
plotted against the tilting parameter y ex = -n & ex /&o in Fig. 4. 
The configuration y ex = -n (i.e. & ex = $0) corresponds to the 
symmetric SQUID potential. The solution of the Schrodinger 
equation for this potential and this coupled SQUID configura- 
tion has two closely spaced energy levels (at the mentioned pa- 
rameter range above) comprising the best pseudo-degenerate 
approximation we ever found to the manifest double degen- 
eracy in Fig. lb. At this point we read from the vertical axis 
in Fig. 4 that Ci = 10 -4 eV. As the double-well potentials 
are tilted in parallel and simultaneously in both SQUIDs the 
shift obtained for each eigenenergy is approximately linear in 
the bias flux & ex . This indicates that the coupling can be 
primarily represented as a first order perturbation in <& ex . A 
comparison of the Fig. 1 with the calculations shown in Fig. 4 
indicates that the coupled rf-SQUID model is in essence an 
isolated four-level quantum system. Recently, Schrodinger 
Cat states have been obtained experimentalUj by manipulat- 
ing these two couplings in the MQT regimefl For the physi- 
cally relevant range < y cx < yt x the centers of the splitting 
energies in Fig. 4 also respect the manifest conservation prin- 
ciple of the energy centers in Fig. 1. In fact, in the numerical 
calculations we found a shift smaller than 2% in the energy 
centers for the whole range < y ex < yt x - 

In this paper we will refer to two different representa- 
tions of the basis states for two coupled SQUIDs. The 
first, which is the even and the odd states of the mutu- 
ally noninteracting SQUIDs, is represented by the vector 
|v) = {|11), 1 00), 1 10), |01> } as in section A. Here in the state 
\mn) m represents the parity state of the SQUID A and n rep- 
resents that of SQUID B. We will designate |v) as the two-bit 
parity basis. The second, which is the set of flux states for the 
coupled and tilted SQUIDs , i.e. both couplings are nonzero, 
is represented by |V) = {| ft). I U>> I IDA it)} which we 
refer to as the two-bit flux basis. 



C. The quantum gate operations 

According to the Hamiltonian (^), the eigen levels in Fig. lb 
corresponding to Ci 7^ and C2 — are the doubly degenerate 
states |v') = U\ |v). Explicitly these are, (|00> + |11))/V5 
and (| 10) + |01))/V2 for E + \Ci\, and, ( 1 10) - |01>)/\/2 and 
(|00) - |11))/V5 for E - \d\. If C 2 is also turned on as 
in Fig.lc, the non-degenerate energy eigenstates become the 
two-bit flux states defined in Eq. (pi) and written explicitly as 
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The choice of the basis in which the fundamental quantum 
gates are to be operated should be determined by the acces- 
sibility of the input and output states by the measurement 
mechanism. This condition actually means that all informa- 
tion in the relative magnitudes and the phases of these states 
should be measurable. The absolute values of the amplitudes 
of the flux states are directly measurable whereas their rel- 
ative phase factors do not couple to the measurement (see 
below). Since it is crucial to be able to read and write the 
relative phase between the flux states, we choose the two-bit 
parity basis |v) instead of the two-bit flux basis |V) to per- 
form the gate operations. The flux basis is used for the I/O 
operations. Before and after performing the logic gates it is 
therefore needed to switch between the |v) and the |V) bases. 
Notice that the matrices U\ and U2 are two-bit Hadamard 
transformations. In particular, it can be verified that 
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is the one bit 



where 1 is the 2x2 unit matrix. The matrix 
Hadamard transformation. The transformation from the nat- 
ural eigenbasis | V) to the parity basis |v) can therefore be per- 
formed by two successive two-bit Hadamard transformations. 
In this way any phase information of a quantum state encoded 
in |v) basis contributes to the relative amplitudes of the |V) 
basis yielding direct accesibility to the measurement. This 

transformation can be summarized as I V) |v') — ^> |v). 

The effective result of this transformation is that each state 
in the parity basis is finaly assigned a distinct eigen energy as 
shown below. This transformation is equivalent to effectively 
assigning each component \mn) of the |v) basis an energy 
eigenvalue. The logic operations are then performed on the 
|v) basis. We now demonstrate the fundamental single-qubit 
operations, i.e. the phase flip and the Hadamard transforma- 
tion. 

Single qubit operations: 

Although the mechanism is based on four states, it is pos- 
sible to identify two individual flux qubits and perform single 
bit operations on each one independently. In the parity ba- 
sis, I ran) we assign the first bit as the control (qubit no.l) 
and the second bit as the target (qubit no. 2). By definition, 
the single qubit operations are performed on the target in the 
four dimensional state space unconditionally from the state 
of the control bit. Whereas the controlled operations (two- 
qubit operations) are performed when the control bit is in a 
particular, i.e (0 or 1) state. With this convention we now 
start with the single qubit operations. 

I. The phase flip is a special case of the phase evolution 
described by 




|H) 
|01) 

e"^ 1 00) 
e"^ 1 10) 



(7) 



The phase flip is obtained in ((?]) at <j> — ty/2 which implies that 
the state of the second bit is changed as 1 — > 1 and 0^—0 
independent from the state of the first bit. It is shown below 
that this can be obtained as a direct sum of two controlled 
phase flips. 

The phase evolution Z in Eq. (Q) is achieved in the two 
coupled SQUID realization in two steps. We switch the cou- 
pling configuration of the system to C2 = Ci which makes the 
states 1 01} and 1 00) degenerate at the major center Eq. By 
this switching, the remaining (nondegenerate) states, i.e. |11) 
and 1 10) are now at the energies Eq + 2C\ and Eq — 2Ci re- 
spectively. The four state system is then permitted to freely 
evolve in this configuration. Degenerate states do not acquire 
any relative phase with respect to each other. At this first 
stage the phases acquired by the states are described by 



(8) 



where (j> = 2Ci t as the time dependent phase. In the second 
stage we set C2 = —C\ which turns the states |11) and 1 10) 
to become degenerate at the energy Eq and 1 01) and 1 00) 
nondegenerate at energies Eo + 2Ci and Eq — 2Ci respectively. 
An identical phase evolution is then performed on the states 
1 01} and 1 00) . The net transformation is given by 
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which, after factoring out the overall phase e' 2 ^ E ° is equiv- 
alent to Eq. (jjj). The phase flip is obtained at (f> = n/2 corre- 
sponding to ti = 7r/4Ci which can be described as 
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-iB t le iH Cl ,- Cl tl (C 2 ^-Cl) 

- iB ° t i e <Hc i' c i tl l(Ci -vCi) . (10) 



where the action of the operators (matrices) is defined right- 
wards. Notice that, Eq. ( |10[ ) is manifestly the composition of 
two controlled phase flips where these two operations com- 
mute. Subtracting the overall phase, the first square bracket 
is equivalent to [| 11} -» |11), 1 10) -> — 1 10}] . This is the 
controlled phase flip 1 — ► 1 and — > —0 when the control 
bit is in the state 1. Similarly, the second bracket denotes 
the same controlled phase flip when the control bit is in the 
state. Hence, this method proves an efficient realization of all 
controlled operations as well (see the double qubit operations 
below) . 

II. The Hadamard ( 



H 



transformation is a special case 
that can be obtained by rotations and phase flips on the Bloch 
sphere and its realization is very similar to the one discussed 
already in the phase flip. The idea is to perform two con- 
sequtive controlled rotations by coupling the non-degenerate 
states to an external pulse resonant at the energy difference 
4Ci . Basically the two-bit Hadamard gate can be obtained in 
the following steps 
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In the first step [C2 — > Ci] the inner states, i.e. 1 01) , 1 00) , 
are made degenerate at the major center Eq. In the second 
step a Hadamard transformation is performed (up to an over- 
all phase e~" r ' 2 ) on the outer states |11) and 1 10) (for those 
the first bit is 1) when the inner states (the first bit is zero) 
are degernerate. In this term, the states 1 1 1) and 1 10) are ro- 
tated by R x (n/4) sandwiched between two phase evolutions 
at <f> = 2 x 7r/4 indicated by R z (7r/4) . The overall phase accu- 
mulated on the four states is e lE ° * 2 , where t2 is the duration 
of the pulse given by ti = 7t/(4k) with k describing the cou- 
pling strength of the levels to the external pulse [The coupling 
mechanism to the external field is similar to the ion trap ex- 
periments and will not be repeated here]. The third factor 
[C2 — > — Ci] indicates that C2 is now switched in the oppo- 
site direction swapping the inner and the outer states. The 
fourth factor in (Jill) is another Hadamard transformation in 
the interchanged configuration (when the first bit is 1). The 
same front factor e lE ° * 2 as in the step 1 also appears here as 
an overall phase. After the fourth step the inner and outer 
states should be swapped back to the initial configuration. 
Note that in all single and double qubit operations there is an 
overall phase acquired by all qubit states. Also note that the 
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gate described by Eq. ([jl]) is a product of the two trans- 
formations [indicated by the square brackets therein] where 
each individual bracket is a controlled H conditioned by the 
state of the first bit. 

Double qubit operations: The controlled phase flip was con- 
sidered previously. It basically comprises the first or the sec- 
ond half of the Z ^ operation in Eq. (p"(i|). We will therefore 
not consider it here again. The other crucial double qubit 
operation is the Controlled-NOT gate. 
Controlled NOT 
Ideal! 
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It is much easier to obtain CNOT in this proposed mech 



anism, which can be done by a single controlled rotation by 
7r/2. This can be written as 



D. Operational time scales and dephasing 

The time scale required for the gate operations can be 
rougly considered to be independent from the particular gate 
operation. We consider the operational time as t op — (AE)^ 1 
where AE describes the energy difference of the states at the 
configuration where the gate operation is performed. For in- 



stance, a typical scale is the phase flip t a 



ti 
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be read off from Fig. 4 at the degeneracy point y ex = y% x . 
We find that Ci ~ 5 10~ 3 eU which yields t x ~ 5ps. The 
switching tjkne is possible for less than Ins between two gate 
operations.^! By comparing these two time scales, we estimate 
that the gate operations are dictated by the switching time. 

As the dephasing mechanisms are concerned certain ex- 
pected sources and their decaherence time scales have been 
previously analyzed in detailBEi SQUIDs in the flux regime 
are known to be insensitive to the charge fluctuations. The 
decoherence time for the charge fluctuations in the flux dom- 
inated regime has been estimatedu to be t v ~ 0.1s. Among 
the other sources there are, the quasiparticle tunnelm* (t v ~ 
lms) and electromagnetic radiation of the junctiortfj (t v ~ 
10 3 s). More substantial dephasing effects are known to arise 
from the near distance dipole-dipole interactions between the 
SQUIDs (t v ~ 0.1ms). 

The relaxation time in the time dependence of the basis 
(flux or parity) states is the most crucial element for the I/O 
operations based on the density matrix measurements. For 
the parameters considered here it can be calculated by the 
explicit expression given by Leggett et al. and Weiss Ref. [11] 
for which we find t re j ~ lfis. 



E. Read-in and read-out 

The state of the coupled squids is observed before and after 
the computation by the flux measurement made by two dc- 
squids each independently facing one SQUID loop as shown 
in Fig. 2. Before any read-in or read-out made, a back trans- 
formation is needed between the operational |v) basis to the 
physical flux basis |V). Suppose that a particular logic opera- 
tion is denoted by the matrix T in the basis |v). The formula 
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CNOT = e 



Ci] Rx(n/2) [C2 — * C2] ■ (13) transformed by, 



describes the effect of T in the basis |V). Here T is any 
desired composition of the single and double qubit operations. 
Suppose T = \z\ k . Using @ we find that the flux states are 



If one ignores the overall phase, Eq. ( |13[ ) amounts to setting 
the coupling so that C2 = Ci, applying the rf-field to induce 
the rotation R x (n/2), and finally switching back to the initial 
configuration of the Ci and C2. The full 
the |v) basis is 
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We give a second example from T ■■ 
to this case we obtain 
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Applying @ 



which amounts to the flip of the second bit |0) 
1) — > i|0) when the first bit is in the 1 state. 



(14) 
and 



|V> 



(1-i) (1 + i) 

1 I (1 + i) 

2 I ~(l + i) (1-i) 
-(1-i) (1 + i) 



IV) . 
(17) 



5 



The unique correspondence between the flux states and the 
parity basis permits read-in and read-out operations to be 
performed by measuring the flux as shown in Fig. 2. The 
flux measurements simulate single event quantum mechan- 
ics. Suppose an output state to be measured is \ip) = di | tT 
)+a2 | it) +^3 | Tl)+ a 4 I ||) A single flux measurement yields 
one of the flux states in the output with the probability \a,i\ 2 
associated with that flux state. The partial amplitudes of the 
flux states can therefore be inferred only after repeated mea- 
surements. In this context, the final measurement is similar 
to the quantum computation model using ion traps. 



F. Discussion 

We demonstrated a new theoretical mechanism for quan- 
tum logic gate operations by manipulating degeneracy of a 
symmetric four level system by two variable coupling param- 
eters. A realization of this system using rf-SQUIDs in the 
flux regime is made in which only one of the couplings is vari- 
able. Physical examples in which both couplings vary are 
very likely to exist but currently unknown to us. In that case 
an ideal realization of the proposed mechanism is possible 
with the particularly added feature of Fig. la namely storing 
the computed information by suppressing the energetic inter- 
ference between the states. The realization of the gates for 
the fundamental single and double qubit operations is also 
demonstrated. As the scalability of this SQUID realization 
is concerned we have less to say at this moment. Masking 
and evaporation techniques that are already well established 
can be used to fabricate arrays of the qubit pairs. The cou- 
pling between different pairs can be realized in a similar way 
as those that are coupled as members of a single pair. We 
continue our current effort in the direction of simulating two 
and three such pairs coupled inductively. 
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FIGURES 

Fig.l: Lifting the four- fold degeneracy stepwise in the pro- 
posed mechanism. In the ordering of the flux states shown in 
(c) [see Eq. (^) for the definition of the flux states] with respect 
to the energy eigenvalues we considered < C% < C\. Note 
that the initial four-fold degeneracy causes the energy centers 
indicated by the horizontal dashed lines to be independent of 
the couplings. 

Fig. 2: Inductive coupling of the two rf-SQUID loops. The 
dc-SQUIDs have two-junctions for balanced read-in and read- 
out the flux at each rf-SQUID loop independently. Although 
the details of the mutual field screening is not shown, the 
circuitry should be designed such that the I/O devices are 
decoupled from each other completely. 

Fig. 3: Pseudo-degenerate energy eigenfunctions and the 
corresponding single well localized flux states. The flux states 
Oil are properly normalized by 1/ \p2. The parameters used 
are fl L = 1.1, y ex =n,C = OApF, L = WOpH. 

Fig. 4: Lifting the degeneracy by tilting the symmetric 
double- well potential in the coupled SQUID pair. The vertial 
scale is the numerically calculated eigen energies in eV units. 
The horizontal variable corresponds to no tilt when y ex = tt 
and critical tilt when y ex = y* x . The tilting parameter y ex /Tr 
simulates the coupling C2 in the proposed mechanism in sec- 
tion A. The parameters used are J3l = 1.1, C = OApF, 
L = WOpH and M = 0.03L. 

Fig. 5: The transformation from the flux basis |V) to the 
parity basis by two Hadamard transformations. After U\ is 
applied on the flux basis, the states are v[ — ( ||) — | || 
))/V2;4 = (|TI) + |IT))A/2; 

va' = (I TT) - I ll))/V2; v'i = (| TT) + I U»/V5 with the 
energies E x = E - Ci - C 2 , E 2 = E - & + C 2 , 
E 3 = E Q +C1-C2, E A = Eq + C1+C2 . 
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